Phaseless three-dimensional optical nano-imaging 

Alexander A. Govyadinovj^] George Y. PanasyukJ^] and John C. Schotlancj^] 

Department of Bioengineering , University of Pennsylvania, Philadelphia, Pennsylvania 

(Dated: June 30, 2009) 

Abstract 

We propose a method for optical nano-imaging in which the structure of a three-dimensional in- 
homogeneous medium may be recovered from far-field power measurements. Neither phase control 
of the illuminating field nor phase measurements of the scattered field are necessary. The method 
is based on the solution to the inverse scattering problem for a system consisting of a weakly- 
scattering dielectric sample and a strongly-scattering nano-particle tip. Numerical simulations are 
used to illustrate the results. 
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The development of tools for three-dimensional imaging of nanostructures is of consid- 
erable current interest [H 121 El Hj- There are multiple potential applications, including 
inspection of semiconductor devices, detection of atoms buried beneath surfaces and char- 
acterization of biologically important supramolecular assemblies, among others. Optical 
methods, especially near-field scanning optical microscopy (NSOM) and its variants, hold 
great promise for nanoscale imaging due to their subwavelength resolution, spectroscopic 
sensitivity to chemical composition and nondestructive nature [5]. Although traditionally 
viewed as a technique for imaging surfaces, near-field microscopy has recently demonstrated 
the capacity to detect subsurface structure [U [3J. Experiments in which a near-field probe 
is scanned over a three-dimensional volume outside the sample suggest that information on 
the three-dimensional structure of the sample is encoded in the data. That is, the mea- 
sured intensity viewed as a function of height above the sample is seen to depend upon the 
depth of subsurface features. However, the intensity images obtained in this manner are not 
tomographic, nor are they quantitatively related to the optical properties of the medium. 

The above noted difficulties have led to the use of inverse-scattering theory to elucidate 
the precise manner in which three-dimensional subwavelength structure is encoded in the 
optical near-field [SI [Tl EJ [91 HOl HH [121 [J3l HH HSl HS1 HZ]- Results in this direction have 
been reported for two-dimensional reconstructions of thin samples [10] and also for three- 
dimensional inhomogeneous media [18]. In either case, solution of the inverse scattering 
problem generally requires measurements of the optical phase, in the form of a near-field 
hologram, an experiment that is notorious for its difficulty The replacement of phase mea- 
surements by phase control of the illuminating field has also been proposed [9]. In this 
approach, the power extinguished from an incident evanescent wave field is used to recon- 
struct the imaginary (absorptive) part of the dielectric susceptibility, leaving the real part 
unrecoverable. 

In this Letter we propose a method for nano-scale optical tomography that relies neither 
on phase-measurements of the scattered field nor on phase-control of the illuminating field. 
Our approach enables the reconstruction of the complex-valued dielectric susceptibility with 
subwavelength resolution in three dimensions. As a proxy for the optical phase, we introduce 
a controlled scatterer, such as an atomic force microscopy (AFM) tip, into the near-field of 
the sample. The power extinguished from the incident field, which illuminates the sample 
and the tip, is then measured in the far-field. Since the tip is placed externally to the sample, 
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changing its position controls the pattern of illumination, which thus modifies the power 
extinguished from the incident beam. The crucial difference from conventional holographic 
techniques is that the interference pattern is regulated by an internal degree of freedom of 
the system (regarded as the sample plus the tip) rather than by external illumination. The 
burden of phase-resolved measurements or phase-controlled illumination is thus replaced by 
the problem of controlling the position of the tip. The readily available nanometer precision 
in probe positioning that is achievable in AFM, in combination with the simplicity of far- 
field measurements of the extinguished power, is expected to allow the practical realization 
of the proposed method. 

We begin by considering an experiment in which a sample is deposited on a planar 
substrate. The lower half-space z < (the substrate) is taken to have a constant index 
of refraction n. The sample occupies the upper half-space z > and is assumed to be 
nonmagnetic. The index of refraction in the upper half-space varies within the sample, 
but otherwise has a value of unity. The upper half-space also contains the tip which is 
placed in the near-field of the sample. The sample and tip are illuminated from below by 
a monochromatic evanescent plane wave and the power extinguished from the illuminating 
field is monitored, as shown in Fig. [T} 

The electric field E in the upper half-space obeys the reduced wave equation 

V x V x E(r) - A; 2 E(r) = 4nk 2 (r,(r) + X (r)) E(r) , (1) 

where r\ is the generally-complex dielectric susceptibility of the sample, x is the susceptibility 
of the tip, ko = 2ir/X is the free-space wave number and the field obeys the necessary interface 
and boundary conditions. The field is taken to consist of two parts, E = Ej + E s , where 
Ej is the incident field and E s is the scattered field. The incident field obeys Eq. in the 
absence of the sample and the tip. The scattered field obeys the integral equation 

E.(r) = k 2 J G(r, r') • E(r') (r,(r') + X (r')) dV , (2) 

where G is the half-space Green's tensor. The power P e extinguished from the illuminating 
field can be obtained using the generalized optical theorem [19J: 

P e = ^Im J E*(r) ■ E(r) fo(p) + x (r)) d 3 r , (3) 

where the integration is performed over the volume V, which contains both the sample and 
the tip. 
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Suppose that the tip is a strongly-scattering, possibly metallic, nanoparticle and the sam- 
ple is a weakly-scattering dielectric. We may then compute the electric field perturbatively, 
accounting for all orders of scattering from the tip and one order of scattering from the 
sample. We find that the resulting perturbation series can be resummed and, neglecting 
contributions arising solely from the sample or the tip, consists of a sum of three terms [IT] . 
The first, or 'TS', term corresponds to scattering from the tip and then from the sample. 
The second, or 'ST', term is due to scattering from the sample and then from the tip. The 
third, or 'TST' term, arises from scattering first from the tip, then from the sample and 
finally from the tip. Note that two additional terms originating solely from scattering by 
the sample or the tip contain no structural information and will be omitted. In practice 
they can be removed by calibration. To proceed further, we must specify a model for the 
tip. We treat the tip as a small scatterer with susceptibility x(r) = a <5( r ~~ r t), where r t 
is the tip's position and a® is its polarizability. Resummation of the perturbation series, 
as explained above, leads to a renormalization of the polarizability of the tip of the form 
a = ao/(l — 2i/c 3 «o/3), which includes the lowest order radiative corrections to the bare 
polarizability but neglects the dependence on the tip height above the interface [I7J 120] . 

It follows from Eq. ^ that the extinguished power can be expressed as a sum of contri- 
butions of ST, TS and TST types: 

P e (r t ) = ^ I Y, KiP) ( r t,r)r ] ^(r)d 3 r , (4) 
p=i 

where r]^(r) = r]^*(v) = r)(r), the kernels K^'(r t ,r) are defined by 

KW(t u t) = -K^\v u v) = aE*(r t ) • G(r t ,r) ■ E,(r) + «E*(r) ■ G(r,r t ) ■ E^) 

+aVE*(r t ) ■ (G(r t , r)G(r, r t )) ■ E<(r t ) (5) 

and the dependence of the extinguished power on the tip position has been made explicit. 

We will assume that the sample occupies the region < z < L and that it is illuminated 
by a plane wave of the form Ej(r) = E exp(iqj • p + k z z). Here r = (p,z) and the field 
has amplitude E , transverse wave vector q^, and k z = a/ (nk ) 2 — qf. The extinguished 
power is measured for a discrete set of tip positions located on a three-dimensional cartesian 
grid with transverse spacing h and longitudinal spacing Az. Note that the tip occupies the 
region L < z < L t and thus does not overlap the sample. 
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It will prove useful to perform a two-dimensional lattice Fourier transform of the sampled 
extinguished power in the plane z = z t , namely P e (q,z t ) = ^ p exp(zq • p)P e (p, z t ). Here 
the sum is carried out over all lattice vectors and q is restricted to the first Brillouin zone 
(FBZ) of the lattice. Next, we require the plane-wave decomposition of the tensor Green's 
function 

G(r, r) = J expfiq ■ (p - f/)]gq(z, z) , (6) 

where the form of g q is given in Ref. [21j . Making use of this result and carrying out the 
lattice Fourier transform, we find that Eq. (|4]) becomes 

P e (q, z t )= R(P) ^ z <> z )^ W (q, *)dz , (7) 

Jo P =i 

where ifW(q; Zt,z) is defined as 



tf«(q; ^ z) = aT (7E* ■ g qi _ q (^, z) ■ E + 7 *E* ■ g ni+cl (z, z t ) ■ E 

f d 2 q' \ 

and K (2) (q; z t , z) = i^ (1) *(-q; z t , z). Here T = ckl/(Aih 2 ) exp[— 2lmk z Zt[, 7 = exjp[ik z (z— z t )) 
and ^(q, z) = J d 2 p exp(iq- p)r)(p, z). Note that for fixed q, Eq. ([7| defines a one-dimensional 
integral equation for fj^ (q, z). 

The inverse scattering problem we consider consists of recovering r}^ p \ for p = 1,2, from 
measurements of P e . This corresponds to solving the integral equation Eq. ([7]). If it is 
known, apriori, that the susceptibility rj is purely real or imaginary, then the inverse prob- 
lem is formally determined and the solution to Eq. ([7]) is readily obtained by singular value 
decomposition (SVD) [9]. However, if rj is complex- valued, the inverse problem is under- 
determined. To resolve this difficulty, it is necessary to introduce additional data, which we 
take to consist of a second set of measurements. That is, two sets of measurements must 
be carried out for each location of the tip, yielding P el and P e2 , one for each incident plane 
wave with transverse wave vectors qi^ and amplitudes E12, respectively. In this manner, it 
is possible to reconstruct both 77 and rj* simultaneously, which is equivalent to recovering the 
real and imaginary parts of 77 from power extinction measurements. Following the approach 
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of Ref. [22], we find that the solution to the integral equation ([7]) is given by the formula 

where i,j = 1,2 label the incident waves. Here is the inverse of the matrix whose 

elements are given by 

M«,-(q; z t , z' t ) = f Y, *i P) (* z <> <> z ) dz • ( 9 ) 

Jo p 

An inverse Fourier transform is then applied to obtain a transversely bandlimited approxi- 
mation to t](r) with bandwidth 2ix/h. 

To demonstrate the feasibility of the inversion, we have numerically simulated the recon- 
struction of 1]{y) for a collection of point scatterers. The left column in Fig. [2] shows the 
configuration of the scatterers. The tip was modeled as a small sphere of polarizability cto = 
(e — l)/(e + 2)R 3 with radius R = 8x 10~ 2 A and permittivity e = — 1 1.39 + 0. 13z, which corre- 
sponds to silver at a wavelength A = 550 nm. The incident fields were taken to be evanescent 
plane waves with transverse wave vectors qi = (3.15/co/vr, 0) and q2 = (0, 3.25/eo/vt), and vec- 
tor amplitudes E 1 = (-0.521, -0.714, 0.468) and E 2 = (-0.507, -0.714, 0.483), respectively. 
The susceptibility rj was reconstructed on a 40 x 40 x 20 cartesian grid whose transverse 
extent was 0.4A x 0.4A and height in the ^-direction was 0.08A. The forward data were 
calculated from Eq. Q for the positions of the tip center located on the same 40 x 40 
transverse grid with 20 steps of size Az = 0.001A in the z direction, beginning 0.16A from 
substrate. The integral in the kernel ^ was numerically evaluated using a trapezoidal rule 
with 300 points in each direction, spanning six Brillouin zones. We note that computation of 
the matrix inverse Mj" 1 requires regularization, which we carry out by retaining only those 
singular values in the SVD of My which are larger than a prescribed cutoff. 

In Fig. [2] we present reconstructions of the real and imaginary parts of 77. Tomographic 
slices are shown in the planes z = 0.016A and z = 0.068A. It can be seen that the scatterers 
in the top layer (nearest the tip) are better resolved than the scatterers in the deeper layer. 
This is due to the decay of high-frequency evanescent waves with depth and is a typical 
feature of tomographic reconstructions in the near-field [7J. It may also be observed that 
the reconstructions of the imaginary part of the susceptibility are of higher quality than 
those of the real part. This effect may be explained by noting that the extinction of power 
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due to absorption is greater than that due to elastic scattering in the near-field, since the 
optical phase changes minimally in the near-zone of the scatterer. 

In conclusion, we have shown that the three-dimensional subwavelength structure of an 
inhomogeneous scattering medium can be recovered from far-field measurements of the 
extinguished power. Remarkably, neither phase control of the illuminating field nor phase 
measurements of the scattered field are required. Our approach is based on the solution 
to the inverse scattering problem for a system consisting of a weakly-scattering sample and 
a strongly scattering nano-scale tip. In principle, the theory can be extended to treat the 
case of a strongly-scattering sample by inversion of an appropriately resummed perturbation 
series, taken to all orders of scattering in the sample pj)]. It is important to note that the 
observed subwavelength resolution is due to the modification of the near-field of the sample, 
resulting from the presence of the tip. It is not due to the illuminating evanescent wave, as 
is the case in total internal reflection tomography [HJ[TT]. The evanescent wave is introduced 
solely for the experimental convenience of measuring the extinguished power. 
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FIG. 1: Illustrating the experiment. The tip scatters the incident evanescent field and modifies 
the interference pattern in the sample, which has dielectric susceptibility rj(r). The power P 
extinguished from the illuminating field is measured as the tip is scanned on a three-dimensional 
grid in the near-zone of the sample. 
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FIG. 2: (color online) The model (left) and simulated reconstructions (right) of Rer/(r) and Imr/(r). 
The scatterers are distributed in two planes at z = 0.016A (top) and z = 0.068A (bottom). The 
scatterers are separated in-plane by 0.05A and 0.2A in the x- and y-directions, respectively. Each 
image is normalized to its own maximum and any small negative values are not displayed. The 
field of view of each image is 0.4A x 0.4A. 
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